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A fast, simple to use Monte Carlo based event generator is presented 
which is intended to facilitate simulation studies and the development of 
analysis software for the Solenoidal Tracker at RHIC (Relativistic Heavy 
Ion Collider) (STAR) experiment at the Brookhaven National Laboratory 
(BNL). This new event generator provides a fast, convenient means for pro- 
ducing large numbers of uncorrelated A+A collision events which can be used 
q ■ for a variety of applications in STAR, including quality assurance evalua- 

tion of event reconstruction software, determination of detector acceptances 
and tracking efficiencies, physics analysis of event-by-event global variables, 
O ' studies of strange, rare and exotic particle reconstruction, and so on. The 

user may select the number of events, the particle types, the multiplicities, 
><J | the one-body momentum space distributions and the detector acceptance 

ranges. The various algorithms used in the code and its capabilities are ex- 
plained. Additional user information is also discussed. The computer code 
implementation is called MEVSIM. 
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1 Introduction 

When carrying out STAR detector simulations and physics analysis studies one 
often needs many events of a similar type that can either be readily selected or 
tightly controlled by the user to meet the specific needs of the simulation project. 
For example, in evaluating the performance of the event reconstruction software 
it is beneficial to use events which have the same particle multiplicities and dis- 
tributions from event-to-event in order to help identify anomalous changes in the 
tracking performance. Events with flat distributions in rapidity facilitate the deter- 
mination of detector acceptances and tracking efficiencies. Events with enhanced 
numbers of multi-strange particles and other rare particles facilitate the study 
of the corresponding reconstruction and physics analysis software. When using 
standard events generators such as hijing, VENUS, RQMD, vni, etc. the user has 
little control over the simulated particle production other than to select the beam 
energy and colliding nuclei, impact parameter, and perhaps the settings of some 



physics process control switches. Furthermore, most of these codes require signifi- 
cant amounts of cpu-time. For many STAR simulation studies these sophisticated 
models are not really necessary and unless the user is interested in the inherent par- 
ticle correlations which a given model may produce, a simple Monte Carlo particle 
generator will generally suffice. 

In response to this need we have developed a standalone Fortran 77 computer 
code called mevsim. The code is designed to provide a fast, convenient means for 
producing large numbers of uncorrelated A+A collision events which can be used 
for a variety of applications in STAR, including quality assurance evaluation of 
event reconstruction software, determination of detector acceptances and tracking 
efficiencies, physics analysis of event-by-event global variables, studies of strange 
and rare particle reconstruction, and so on. The user may select the detector 
acceptance ranges, the number of events, the particle types, the multiplicities, and 
the one-body momentum space distributions with respect to transverse momentum 
(Pt), rapidity (y) and azimuthal angle (0) from a menu of models. 

The particle multiplicities can either be fixed for each event or allowed to 
randomly vary from event-to-event according to Poisson statistics. The parameters 
of the one-body momentum space distributions (e.g. slope parameter, rapidity 
width, radial expansion velocity) can either be fixed or randomly varied from 
event-to-event according to a Gaussian distribution in order to simulate dynamical 
fluctuations. In addition the user may specify a reaction plane angle and include 
anisotropic flow. The overall multiplicities for all particle types in each event can 
also be scaled by a common, multiplicative factor in order to simulate trigger 
fluctuations. In addition to the usual list of particle types in geant, the code also 
includes a number of short lived resonances, several of which have explicit mass 
distributions. These resonances are not decayed. Finally, the events produced are 
assumed to be in the A+A center-of-momentum frame which for symmetric beam 
species at RHIC coincides with the lab frame. 

In the following sections we describe the random sampling methods used in 
the code, the one-body momentum space distribution models, the acceptance cuts 
for application to STAR, the implementation of resonances, trigger fluctuation 
simulations, reaction plane and anisotropic flow simulations, and the standalone 
code and related files. 

2 Random Sampling Methods 

The method used in the code to randomly sample a one-dimensional probability 
distribution is described first, followed by an explanation of the random selection - 
random rejection algorithm which is used to populate momentum space according 
to the specified one-body distributions. 

In order to randomly sample a one-dimensional probability distribution P(x) 
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throughout the domain X\ < x < £2 one first obtains the integral, 

F(x) = [ X P{x')dx' , (1) 



normalized such that F(x2) = 1. By finding the value of x such that F(x) = f, 
where / is a randomly selected number between and 1, a random collection of 
x values can be obtained. A histogram of these values of x will be statistically 
consistent with the probability distribution P(x). This method is much faster 
than the usual random selection - random rejection method since cpu-time is not 
wasted on trial values of x which end up being discarded. In the code the above is 
implemented by integrating P(x) on a fixed mesh, randomly selecting / and inter- 
polating the inverse F function on an uneven mesh using a Lagrange interpolation 
routine. 

The preceding method is used to sample the reaction plane angle, the overall 
multiplicity rescaling factor, the particle multiplicities, the parameters of the one- 
body momentum space distribution model, the anisotropic flow parameters, and 
the resonance masses. Generally these follow Gaussian probability distributions, 

P(x,x mean ,cT) = iVexp -(1 - x mean ) 2 /2(J 2 , (2) 

where N is a normalization constant and the user specifies the mean value, a and 
number of as (domain) over which x is to be sampled. The particle multiplicities 
follow a Poisson distribution, 

P(x, X me an) ' j (3) 

x\ 

with mean value x mean and standard deviation y/x meajI1 . The user specifies the mean 
multiplicity and the (±) number of standard deviations about the mean over which 
to sample. The resonance masses follow approximate Breit-Wigner distributions 
which are listed in Sec. [5| In order for each of these quantities to fluctuate the user 
must input positive values for both the standard deviation and the (±) number of 
standard deviations over which to sample. 

The momentum space distributions, being two-dimensional in p? and y, are 
sampled by the usual random selection - random rejection method. Furthermore, 
since the anisotropic flow parameters depend on and y, it was more convenient 
to also sample the azimuthal {(f)) distribution by the random selection - random 
rejection method, rather than use the preceding integration-interpolation method. 

The arbitrary, two-dimensional function g(u,v), defined in the domain Ui < 
u < Ui and v \ < v < V2, is randomly sampled as follows. The maximum value of g 
(fi'max) m the domain is determined. Values of u and v in the domain are randomly 
selected and the quantity g(w, t>)/g max is computed. A random number from to 
1 is selected and if it is < g(u,v)/ g max (where g must be non-negative which is 
the case for particle densities) , then the randomly selected coordinates u and v are 
retained; if not, these coordinates are rejected and another trial set of u and v is 
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selected. Such methods are reliable but can suffer from poor cpu performance if 
the function g(u,v) is strongly peaked within its domain. For RHIC collisions the 
Pt and y distributions are not expected to have large spikes and we expect this 
sampling method to perform well. However, as noted in Sec. [7], care should be 
exercised when including px and y dependence in the anisotropic flow in order to 
avoid poor cpu performance. 



3 Momentum Distribution Models 

For all momentum space distribution models used in the code the (pr, y) and 
azimuthal (0) dependences are factored. The various models for the {pr, y) depen- 
dences are described in this section, while the dependence models are explained 
in Sec. [j]. 

A total of six (p T , y) distribution models can be selected; the same model must 
be applied to each particle ID type in a run. The models include the following: 



1. Factorized pt and y dependence 

d 2 N 



dp T dy 



Ap T e- m ^ T e- y2 / 2(72 y (4) 



where A is a nomalization constant, my = ym 2 + p|, m is the particle 
mass, T is the m T slope parameter or "temperature," and a y is the rapidity 
distribution width. 

2. Non-expanding spherical thermal source (Pratt) ||, 

pL = A PT Ee-^ T (5) 

where E = m<rcosh(?/) is the total energy. 

3. Non-expanding spherical thermal source (Bertsch) [[|, 

d 2 N Ap T E 
dp T dy e E / T — 1 ' 

4. Spherically expanding, thermally equilibrated source (Pratt) H, 



(6) 



dp T dy 

x 



sinh(, p ) + T f smh(y p ) _ 



(7) 



where y v = jvp/T, 7 = — v 2 , v is the radial expansion velocity param- 

eter, and p = \/E 2 — m 2 is the total momentum. 
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5. Arbitrary, factorized and y dependence loaded bin-by-bin, 

d 2 N 



Af(pr)9(y) (8) 



dprdy 

where / and g are input histograms. 
6. Arbitrary, two-dimensional function of pt and y loaded bin- by-bin, 

-^ r = Ah(p T ,y) (9) 
dp T dy 

where h is a two-dimensional histogram. 

4 Kinematic Acceptance 

The generated particles are confined within specified kinematic boundaries. The 
acceptance cuts are designed to be applicable to the STAR geometry. The ac- 
ceptances are defined by lower and upper values for p T , pseudorapidity (77) and 
azimuthal angle ((f)) (between 0° and 360°). 

5 Resonances 

The available particle types include those listed in geant || plus a number of 
additional particles and resonances which are listed in Table | along with their 
assumed mass and width parameters. For most of these the widths are sufficiently 
broad to warrant including the mass fluctuations explicitly in the generated particle 
list. The mass distribution models assumed in the code are described in this 
section. Since the mass and width parameters depend somewhat on the model, 
the parameters listed in Table p] may not agree exactly with similar values tabulated 
elsewhere, such as in Ref. f7|. In the code the resonance masses are sampled after 
the pt, y and distributions have been determined. Hence, in the final output 
particle list the p?, y, <t> and mass distributions will be exactly as specified by 
the user. However, because the masses are changed after pt and y are fixed, 
the pseudorapidity values will change and the pseudorapidity distributions will be 
altered by the mass fluctuations. 

The following resonance shapes (un-normalized) were assumed: 

1. p-meson || 

W{M) = (M-A^ + rg/4 • (10) 

where q is the cm. momentum in the final it — it system, and M p and T p are 
listed in Table |I[ 
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2. c<> meson 



M 2 r 

W(M) = — , (11) 

1 ' {M 2 - M 2 ) 2 + M 2 Tl{s) ' 1 ' 



where F u (s) = r w (M/M w ) and and T u are given in Table |. 
3. 0-meson jnj 

"w- w-SS+w iB)'- (12) 



where q is the cm. momentum in the final K — K system, E — y g 2 + M^, 
M K is the kaon mass, q and E are evaluated for M = M^, and and 
are given in Table |. 



4. A-resonance jnj 

2(gfi) 3 r A 

r(g) = TTW (14) 

where gi? = 0.81g/m,r, g is the cm. momentum in the final it — N system, 
and \i is the reduced, total energy of the 7r — N system given by, 

M = , = \I<1 2 + ml,e N = sjq* + M% , (15) 

Cvr + (-N 

where and Mjv are the pion and nucleon masses, respectively, and Ma 
and Ia are given in Table |. 



5. K* resonances [12 



W{M) = (M 2 -M 2 K ,) 2 + T(q) 2 M 2 K , ' (16) 
2iy( g /g ) 3 

r(g) = TT(^F' (17) 

where q is the cm. momentum in the final 7r — K system, go is evaluated for 
M = Mk* i and Mk* and Tk* are given in Table III 



6 Trigger Fluctuations 

In the STAR A+A event stream, for which the Level triggers are based on inte- 
grated yields in the central trigger barrel (CTB) and multi-wire chambers (MWC), 
the actual particle multiplicities would fluctuate from event-to-event, even if the 
trigger acceptance settings could correspond to an ideal 5-function. To simulate 
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this effect in the data stream we included in the code the ability to rescale the 
multiplicities of all particle types in the event by a common multiplicative fac- 
tor, /. This factor is selected randomly for each event according to a Gaussian 
distribution which the user specifies. 

Ideally, the nominal mean multiplicities for each particle type, iVpid, should be 
multiplied by / and the resulting Poisson distribution with mean value fN pi ^ and 

standard deviation J fN pid sampled to yield the actual particle multiplicities for 
the event. However, it turns out to be much more straightforward in the code and 
to have less impact on the cpu demand, to instead, sample the nominal multiplicity 
distributions first and then scale the results to get the actual particle multiplicities 
for the event. Using this method the particle multiplicity distributions will have 
the correct mean value of fN pi d but the incorrect (non- Poisson) standard deviation 
°f f\jN P id which is in error by the factor \ff. It is anticipated that / will be within 
a few percent of 1.0 such that the actual, non- Poisson multiplicity distributions 
will only err by a few percent. 



7 Reaction Plane and Anisotropic Flow 



The detection and study of directed flow, elliptical flow and perhaps higher modes 
of particle flow in non-central A+A collision events at RHIC are of great interest to 
the STAR collaboration. These phenomena can be simulated with this generator. 
The user may select from among four options for the reaction plane: 

1. Ignore the reaction plane, events are azimuthally invariant, i.e. no flow. 

2. Assume a fixed reaction plane angle; useful for testing analysis code and 
debugging. 

3. Assume a randomly fluctuating reaction plane angle which follows a Gaus- 
sian distribution; useful for evaluating reaction plane measurement resolution 
effects. 

4. Assume a random reaction plane angle, uniformly distributed between 0° and 
360°, which is the way the data will be acquired in the experiment. 



The model assumed for the anisotropic flow is from Ref. [13|, where 



E 



d 3 N 



d 2 N 



dp 3 2irpT dprdy 



1 + ]T2K(pt,?/) cos[n(0 - 



n=l 



where ^> r is the reaction plane angle, is the azimuthal angle of the outgoing 
particles, and the parameters V n (pT,y) are given by [II ], 



V n (pT,y) = (V nl + V n2 pl)exp (-V n3 (y - y cm 
V n {PT,y) = {V nl + V n2 p T )sign(y 

?/cm) ( Vn3 + Vn4 | y 2/< 



even 



(19) 

,n = odd(20) 
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where sign(y) is ±1 for y > or y < 0, respectively, and y cm is the cm. rapid- 
ity which is for symmetric A+A collisions. Up to six Fourier components are 
permitted; directed flow (n = 1) and elliptical flow {n = 2) are most often used. 

The user should be aware that careless use of the V n (pT,y) parametrization 
could have negative impact on the code's performance. Large values of V n ^ {i = 
1,2,3,4) and/or large pt and y acceptance ranges can lead to very large, unreal- 
istic anisotropy, particularly near the edges of the acceptance. This can result in 
"spikey" <p dependences that will cause the random rejection process to perform 
inefficiently The V n (pT,y) parametrization should therefore be used with cau- 
tion and the user should verify that the "flow" remains within reasonable bounds 
throughout the (p T , y) acceptance range. 



8 Standalone Version 

The source code and other documentation files are available from the authors 
(ray@physics .utexas . edu). Two example input files for the standalone version of 
the code have been designed to reproduce the major particle distributions predicted 
by VENUS and VNI for 200 A-GeV Au+Au central collisions. The input parameters 
are defined and explained in Appendix A. Although resonances may be included, 
they are not decayed in the standalone version. 

The standalone code produces a log file, mult_gen.log, containing diagnostic 
information and error messages if any. For normal, successful runs this file should 
be empty. An output file, mult _gen. out, is produced which includes a lengthy 
header containing the input parameters plus several computed quantities along 
with the list of particles. All events are accumulated in this one file, each be- 
ing separated by event header lines. This output file format follows the gstar 



text format convention described in Ref. [15]. An additional, optional output file 



mult_gen.kin, can be activated which will write out just the list of particles with 
only their geant particle ID type, mass and kinematic information included. This 
is useful for analyses with PAW since this ASCII text file can be easily loaded into 
PAW-ntuples. 



9 Summary 

A fast, convenient Monte Carlo event generator, mevsim, has been developed for 
STAR simulations. With this code the user has complete control over the particle 
types, the multiplicities and the momentum space distributions of the particles in 
the simulated events. Events generated by mevsim are already in daily use as part 
of the event reconstruction quality assurance program in STAR, mevsim events 
can also be used in many other types of simulation and physics analysis projects 
in STAR for years to come. 
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Table I: Resonances included in mevsim. 



Resonance 


ID 


Mass (GeV) 


Width (GeV) 


P+ 


151 


0.783 


0.177 


P~ 


152 


0.783 


0.177 


P° 


153 


0.783 


0.177 




154 


0.782 


0.0084 


i 


155 


0.95750 


0.000208 





156 


1.0194 


0.00442 


T / / 


157 


3.09693 


r\ AAAAAA 

0.000068 


A" 


158 


1.232 


0.071 


A° 


159 


1.232 


0.071 


A+ 


160 


1.232 


0.071 


A++ 


161 


1.232 


0.071 


K*+ 


162 


0.89183 


0.0498 


K*- 


163 


0.89183 


0.0498 


K*° 


164 


0.89610 


0.0505 



A Appendix 

This appendix contains a description of the input for the standalone version of the 
Monte Carlo event generator, multiplicity_gen.f 

(1) n_events - Selected number of events in run. Can be anything 

■ge. 1. 

(2) n_pid_type - Number of particle ID types to include in the 

particle list. e.g. pi(+) and pi(-) are counted 
separately. The limit is set by parameter npid 
in the accompanying include file 'Parameter_values . inc J 
and is presently set at 30. 

(3) model_type - equals 1,2,3,4,5 or 6 so far. See comments in 

Function dNdpty to see what is calculated. 

The models included are: 
= 1, Factorized mt exponential, Gaussian rapidity model 
= 2, Pratt non-expanding, spherical thermal source model 
= 3, Bertsch non-expanding spherical thermal source model 
= 4, Pratt spherically expanding, thermally equilibrated 
source model. 

= 5, Factorized pt and eta distributions input bin-by-bin. 
= 6, Fully 2D pt,eta distributions input bin-by-bin. 
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NOTE: model_type = 1-4 are functions of (pt,y) 

model_type = 5,6 are functions of (pt,eta) 

(4) reac_plane_cntrl - Can be either 1,2,3 or 4 where: 

= 1 to ignore reaction plane and anisotropic flow, 
all distributions will be azimuthally symm. 

= 2 to use a fixed reaction plane angle for all 
events in the run. 

= 3 to assume a randomly varying reaction plane 
angle for each event as determined by a 
Gaussian distribution. 

= 4 to assume a randomly varying reaction plane 
for each event in the run as determined by 
a uniform distribution from to 360 deg. 

(5) PSIr_mean, PSIr_stdev - Reaction plane angle mean and Gaussian 

std.dev. (both are in degrees) for cases 
with reac_plane_cntrl = 2 (use mean value) 
and 3. Note: these are read in regardless 
of the value of reac_plane_cntrl . 

(6) MultFac_mean, MultFac_stdev - Overall multiplicity scaling factor 

for all PID types; mean and std.dev.; 
for trigger fluctuations event-to-evt . 

(7) pt_cut_min,pt_cut_max - Range of transverse momentum in GeV/c. 

(8) eta_cut_min,eta_cut_max - Pseudorapidity range 

(9) phi_cut_min,phi_cut_max - Azimuthal angular range in degrees. 
(10) n_stdev_mult - Number of standard deviations about the mean value 

of multiplicity to include in the random event-to- 
event selection process. The maximum number of 
steps that can be covered is determined by 
parameter n_mult_max_steps in the accompanying 
include file 'Parameter_values . inc J which is 
presently set at 1000, but the true upper limit for 
this is n_mult_max_steps - 1 = 999. 
smp - Same, except for the "Temperature" parameter. 
Lgma- Same, except for the rapidity width parameter, 
cpvel - Same, except for the expansion velocity parameter, 
sir - Same, except for the reaction plane angle 
l - Same, except for the anisotropy coefficients, Vn. 
lltFac - Same, except for the multiplicity scaling factor. 
;s - Number of mesh points to use in the random model 
parameter selection process. The upper limit is 
set by parameter nmax_integ in the accompanying 
include file 'Parameter_values . inc' which is presently 
set at 100, but the true upper limit for n_integ_pts 



(11) 


n_ 


.stdev 


(12) 


n. 


.stdev 


(13) 


n. 


.stdev 


(14) 


n_ 


.stdev 


(15) 


n_ 


.stdev 


(16) 


n_ 


.stdev 


(17) 


n. 


.integ 
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is nmax_integ - 1 = 99. 

(18) n_scan_pts - Number of mesh points to use to scan the (pt,y) 

dependence of the model distributions looking for 
the maximum value. The 2-D grid has 
n_scan_pts * n_scan_pts points; no limit to size of 
n_scan_pts . 

(19) irand - Starting random number seed. 
FOR MODEL.TYPE = 1,2,3 or 4: 

Input the following 17 lines for each particle type; repeat these 
set of lines n_pid_type times: 

(a) gpid - Geant Particle ID code number 

(b) mult_mean,mult_variance_control - Mean multiplicity and 

variance control where : 
mult_variance_control = for no variance in multiplicity 
mult_variance_control = 1 to allow Poisson distribution for 

particle multiplicities for all events. 
Note that a hard limit exists for the maximum possible 
multiplicity for a given particle type per event. This is 
determined by parameter f actorial_max in accompanying include 
file ' common_f acf ac . inc ' and is presently set at 10000. 

(c) Temp_mean, Temp_stdev - Temperature parameter mean (in GeV) 
and standard deviation (Gaussian distribution assumed) . 

(d) sigma_mean, sigma_stdev - Rapidity distribution width (sigma) 
parameter mean and standard deviation (Gaussian distribution 
assumed) . 

(e) expvel_mean, expvel_stdev - S. Pratt expansion velocity 
(in units of c) mean and standard deviation (Gaussian 
distribution assumed) . 

(f) Vn_mean(k) ; k=l,4 - Anisotropic flow parameters, mean values 

for Fourier component n=l. 

(g) Vn_stdev(k) ; k=l,4 - Anisotropic flow parameters, std.dev. 

values for Fourier component n=l . 

Repeat the last two lines of input for remaining Fourier 
components n=2,3...6. Include all 6 sets of parameters 
even if these are not used by the model for Vn(pt,y) (set 
unused parameter means and std.dev. to 0.0). List 4 values 
on every line, even though for n=even the 4th quantity is 
not used. 
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FOR MODEL_TYPE = 5 input the following set of lines for each particle 

type; repeat these n_pid_type times. 

(a) gpid - Geant Particle ID code number 

(b) mult_mean,mult_variance_control - Mean multiplicity and 

variance control where : 
mult_variance_control = for no variance in multiplicity 
mult_variance_control = 1 to allow Poisson distribution for 

particle multiplicities for all events. 

(c) pt_start, eta_start - minimum starting values for pt, eta 

input for the bin-by-bin distributions. 

(d) n_pt_bins, n_eta_bins - # input pt and eta bins. 

(e) delta_pt, pt_bin - pt bin size and function value, repeat for 

each pt bin. 

(f) delta_eta, eta_bin - eta bin size and function value, repeat 

for each eta bin. 

(g) Vn_mean(k) ; k=l,4 - Anisotropic flow parameters, mean values 

for Fourier component n=l. 

(h) Vn_stdev(k) ; k=l,4 - Anisotropic flow parameters, std.dev. 

values for Fourier component n=l . 

Repeat the last two lines of input for remaining Fourier 
components n=2,3...6. Include all 6 sets of parameters 
even if these are not used by the model for Vn(pt,y) (set 
unused parameter means and std.dev. to 0.0). List 4 values 
on every line, even though for n=even the 4th quantity is 
not used. 

NOTE: The pt, eta ranges must fully include the requested ranges 
in input #4 and 5 above; else the code execution will stop. 

Also, variable bin sizes are permitted for the input distributions. 

Also, this input distribution is used for all events in the run; 
no fluctuations in this "parent" distribution are allowed from 
event-to-event . 

FOR M0DEL_TYPE = 6 input the following set of lines for each particle 

type; repeat these n_pid_type times. 

(a) gpid - Geant Particle ID code number 
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(b) mult_mean,mult_variance_control - Mean multiplicity and 

variance control where : 
mult_variance_control = for no variance in multiplicity 
mult_variance_control = 1 to allow Poisson distribution for 

particle multiplicities for all events. 

(c) pt_start, eta_start - minimum starting values for pt, eta 

input for the bin-by-bin distributions. 

(d) n_pt_bins, n_eta_bins - # input pt and eta bins. 

(e) delta_pt - pt bin size, repeat for each pt bin. 

(f) delta_eta - eta bin size, repeat for each eta bin. 

(g) i, j ,pt_eta_bin(i, j) - read pt (index = i) and eta (index = j) 

bin numbers and bin value for full 2D space. 

(h) Vn_mean(k) ; k=l,4 - Anisotropic flow parameters, mean values 

for Fourier component n=l. 

(i) Vn_stdev(k) ; k=l,4 - Anisotropic flow parameters, std.dev. 

values for Fourier component n=l . 

Repeat the last two lines of input for remaining Fourier 
components n=2,3...6. Include all 6 sets of parameters 
even if these are not used by the model for Vn(pt,y) (set 
unused parameter means and std.dev. to 0.0). List 4 values 
on every line, even though for n=even the 4th quantity is 
not used. 

NOTE: The pt, eta ranges must fully include the requested ranges 
in input #4 and 5 above; else the code execution will stop. 

Also, variable bin sizes are permitted for the input distributions. 

Also, this input distribution is used for all events in the run; 
no fluctuations in this "parent" distribution are allowed from 
event-to-event . 
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